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^ il Aging refers to the property of two-time correlation functions to decay very slowly 

Q I on (at least) two time scales. This phenomenon has gained recent attention due to 

^ , experimental observations of the history dependent relaxation behavior in amorphous 

' c/3 ' materials ( "Glasses" ) which pose a challenge to theorist. Aging signals the breaking of 

>.: 

"^ ' time-translational invariance and the violation of the fluctuation dissipation theorem 



> 



X 



during the relaxation process. But while the origin of aging in disordered media 
is profound, and the discussion is clad in the language of a well-developed theory, 



f^ ' systems as simple as a random walk near a wall can exhibit aging. Such a simple 

(N ' 

^D ' walk serves well to illustrate the phenomenon and some of the physics behind it. 
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I. INTRODUCTION 

Aging is an important phenomenon observed experimentally in glassy materials where 
relaxation behavior depends on the previous history of the system |]l|. As a model for a glass 
Edwards and Anderson [Q] introduced an Ising system where the uniform coupling J > 
between neighboring spins is replace by random numbers Jij, drawn from some distribution, 
which are placed on each bond. The properties of such a disordered system are quite different 



when the couphngs are drawn from, say, a gaussian distribution with zero mean, compared 
to those of the ferro-magnetic Ising model with uniform couphngs J > 0: While (for d >2) 
the uniform Ising model reaches an ordered phase below a critical temperature Tc with 
a nonzero spontaneous magnetization as order parameter, no such order emerges in the 
random bond model. Yet, below a certain temperature Tg, the "glass transition," more and 
more spins loose their mobility and freeze into place, but without any collectively preferred 
direction. Thus, while no macroscopic order parameter emerges, microscopically the state 
of the system is ultimately highly (auto-) correlated. The relaxation process towards such 
a state is naturally extremely slow due to the inherent frustration created by the bond 
distribution, and each spin has to "negotiate" its orientation with ever distant neighbors to 
further lower the collective energy within their domain. 

Consider such a spin glass quenched below the glass transition Tg at time ( "age" ) t^ = 
in the presence of a magnetic field. Throughout the sample, domains of various pure states 
develop that grow with characteristic size L(ta) 0. After waiting a time ta = t^, the mag- 
netic field is turned off and the system's response in the form of its remnant magnetization 
is measured. Initially, the response function is only sensitive to the pure, quasi-equilibrium 
states in their distinct domains. But after an additional time scale related to the waiting 
time, tu,, the response spans entire domains and slows down when it experiences the off- 
equilibrium state the sample possesses as a whole. In either regime, though, the remnant 
magnetization decays very slowly, similar to the sequence of graphs (see Ref. |l|]) in Fig. 2 
below. It has been suggested 0| that a simple scaling form might hold for glasses; e.g. that 
the autocorrelation function for the remnant magnetization at the age t^ = t + t^ of the 
experiment is given by 

C{t + t^,t^)=t-^fit/t^) , (1.1) 

where the scaling function / is constant for small argument, and falls either like a power 
law (indication of a single aging time scale) or like a stretched exponential (indication of 
multiple aging time scales [^) for large argument. In any case, time-translational invariance 



in this two-time correlation function C is broken due to the fact that the relaxation process 
becomes dependent on its previous history in form of the waiting time t^. 

Similar quantities have been observed in a variety of phenomenological ||^ and theoretical 
|7HT0| model systems as well, and a debate is raging on the ingredients that are required 



for a model to capture the salient features of glassy materials |rT|. While (as we will see) 



it is easy to produce aging behavior, the more intricate results from temperature cycling 
experiments [l^ are much harder to reproduce, and appear to be a much more stringent 



indicator of the complicated phase-space structure (dubbed "rugged landscapes" [^) real 
glasses possess. 

In fact, a simple random walk near a (reflecting or absorbing) wall is well suited to 
describe the domain-growth picture used above to describe the observed aging behavior. 
Consider a walker starting near the wall n = at t^ = 0. After a waiting time ta = tu, 
she has explored a domain of size Lit^,) ~ t^/^ off the wall. The walk now represents the 
correlation function for a spin at site n in the above spin glass model, and the wall is the 
edge of the (albeit one-sided) domain. In place of magnetization, we measure the probability 
P{n^ t + t^|n, tyj) for the walker to return to the same site n she found herself on at time t^. 
Without the wall, P is of course invariant to shifts in space and time, i. e. t^ is irrelevant and 
no aging behavior can be expected. In the presense of the wall, the walker will venture from 
the site n for small times t after tyj and find herself unconstraint (the "quasi-equilibrium" 
state), but when t ^ t^ she is likely to encounter the wall and carry a memory of that 
constraint back to the site n for t ^ t^. In the next section, we will solve this model 
for a suitably defined two-time correlation function which indeed shows the expected aging 
behavior. 

Of course, in this model it is the wall that explicitly breaks the symmetries of the system 
instead of the dynamics of the process, and the observed aging behavior appears trivial. 
On the other hand, the dynamics of the stochastic annihilation process A + A —>■ 0, which 
is similar to a walk at a wall, yields identical results to those reported below il^, and 



the distinction between explicit and dynamical symmetry breaking is not so clear anymore. 



Finally, the case of an absorbing wall is particularly revealing and illustrates the the effect 
of a process that is slowly dying out (i. e. its norm decays), while true aging should only be 
associated with intrinsic properties of a process that will sustain itself. This point has lead 
to some confusion recently. 

II. THE RANDOM WALK MODEL 

In this section, we will calculate the conditional probability P(n, t|r2o, to) for a walker 
to reach a site n at time t, given that she was at site uq at some previous time to < t in 
the presense of either an absorbing or a reflecting wall. Then, we will compute a simple 
two-time correlation function (see e. g. Ref. 0] for a similar definition) 

C{t + t^, t^) = Y. Pin, t + t^\n, t^)P{n, t^|0, 0) (2.1) 

n 

for a walker to return to a site at time ta = t + 1^, given that she was at the same site at 
time tyj after the start of the walk at time t^ = near the wall ra = 0. For both boundary 
conditions, we find that the walk ages, i. e. shows a scaling behavior according to Eq. ( p,.!] ). 
Note that the breaking of spatial invariance in P due to the wall is crucial for the breaking of 
time-translational invariance in C: For an unconstraint walk P would be invariant in space 
and time, and we would find P{n,t + t^|n, t^^,) = P(0, t|0, 0), and with I]n-P("-)^w|0) 0) — 1; 
it is C{t + tw,tw) = C(t, 0) = P(0,t|0,0), independent of t^. In the presense of the wall, 
spatial invariance is broken while time invariance for P still holds, and C merely simplifies 
to C{t + tw, tw) = J2n -P(^) t|ra, 0)P{n, tw\0, 0) which remains a function of t^. 

To simplify the algebra, we consider instead of the walk equation for P the potential 
problem 

dt(Pir,t) = d^Mr,t), (r>0,t>0), 
0(r,O) = 5{r-ro), 
0(O,t)=O or 9^0(O,t) = O, (2.2) 



where the two boundary conditions in the last hne correspond to an absorbing and a reflecting 
wall, respectively [0. We solve for 0, and identify P{n,t\rQ,0) = (pro{r,t) and 



/■OO 

C{t + tyj, tyj) = dr (j)r{r, t) (f)o{r, t^). 
Jo 



(2.3) 



The Eq. ( [^.2D is easy to solve by converting to an ordinary differential equation (ODE) 
in r using a Laplace transform in t: 0(r, s) = J^ dte~'^^ (f){r,s). The ODE has simple 
exponential solutions in two regions, < r < tq and r > tq, whose four unknown constants 
are determined by the two boundary conditions at r = and r = oo, and by the two 
matching conditions at r = tq where (p is merely continuous. We find 

~ , . 1 



^ro[r,s) 



v^ 



g-v^k-nil _|_ g-v^(»-+T-o) 



(2.4) 



which is easily inverted using standard tables for Laplace transforms [|T^ 

1 



(f>roir,t) 



Tit I 



e 4t ± e 4t 



(2.5) 



In each case, the upper sign refers to refiecting boundary conditions, and the lower sign 
refers to the absorbing case. 

A. Reflecting Boundary Conditions 



Here we insert the appropriate forms of (f) in Eq. ( p.5|) , using the upper sign case, into 
Eq. ( p.3|) and choose rg = for convenience [p!7| : 

2 r<xi 



L^ [Z ~r T^W} ^n 



T^Vttu 



dr { I + e t I e 4*™ 



^ t'h(^] with /(a;) = l + 



TT 



X 



4: + X 



(2.6) 



Thus, while the two-time correlation function does show the aging behavior according to 
Eq. (|1.1|) , its scaling function is particularly trivial with f{x ^ 1) ~ 1 and f{x ^ 1) ~ 2, 
see Fig. 1. 



B. Absorbing Boundary Conditions 



Again, we insert the appropriate forms of (j) in Eq. (|2.5| ), using the lower sign case, into 
Eq. (|2.3|) . But with absorbing boundary conditions, putting the starting position at tq = 
would be instantly fatal for the walker. (In simulations we use n = 1 as starting point.) 
Instead, we can choose tq arbitrarily small and expand to leading order. Since the starting 



point tq is irrelevant for the asymptotic behavior (at large times) considered here |]T7 
can be sure that higher-order corrections in tq will have to be subdominant: 



we 



O l^t + t-ui) ''li 



1 



TTy/tt-u, Jo 



dr 
dr 



1-e-- 



(r-ro) 



4*11, — e 4'™ 



ro 2 



e t 



re ''*™ 



67,! 71 \t^in 



with f{x) 



i + f 



(2.7) 



In this case, we find more interesting scaling behavior with /(a; <^ 1) ~ 1 and /(a; ^ 1) ~ 
4/x. 

But this observed aging behavior does not consider the effect that a walk actually dis- 
appears when reaching the wall which diminishes the norm of the distribution (0 or P). 
Rather, to obtain the intrinsic properties of an infinite walk near an absorbing wall, we 
have to properly normalize the correlation function. To that end, we consider the two-time 
correlation function C{t + t^,t^|6') for a walk that reaches the wall (and disappears) ex- 
actly at time ta = 0, and its generic relation to the intrinsic two-time correlation function 

(0 ie<t^ + t), 

C{t + t^,t^\e) = l (2.8) 

These quantities are related to the two-time correlation function given in Eq. (|2.7| ) (which 
is usually the one that is simulated by averaging over walks of any length up to some cut-off 
in time): Given the probability Pt{6) ~ d^^, t > 1, for the walker to reach the wall for the 
first time at t^ = ^ (at which point the walk disappears without further contributing to the 
statistics in the numerical simulation), we have the identity 
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POO 

cit + t^,t^)= dec{t + t^MO)Pt{e) 

Jo 

--C''''\t + t^,t^){t^ + tf-\ (2.9) 



and thus 

c"" ''(t + tu,, tu,) ~ c{t + tu,, tu,) (tu, + ty 



1 . / 1\ /^ t^^-^ 



r\^ 



Hence, the correct scahng function for the aging behavior of the intrinsic process is given by 

r'^ix)^fix)ii + xy-\ (2.11) 

with /™*''(x <^ 1) ~ 1 and p^^^{x ^ 1) ~ x"^^^. Of course, r = 3/2 from the famihar 
first-passage time of a random walk [|18|, and aging remains intact ahhough the cross-over 
in Z™*'' is less dramatic then before for /. In Fig. 2 we plot results for (7™*"^ from numerical 
simulations, and in Fig. 3 we plot the scaling /™*'' for the data in Fig. 2. 

As mentioned before, the stochastic annihilation process A + A ^ is closely related to 
the random walk model with an absorbing wall and, indeed, the intrinsic scaling behavior 
found here (aside from an overall factor of v^) coincides with the on reported in Eq. (9) of 
Ref. [H (for t ^ t + t^ and ^ -^ t/{t + T^)). 

While in this walk model the aging behavior remains intact even for the intrinsic prop- 
erties of the process, it is important to note that in some cases the observed aging behavior 
can be entirely attributed to improper normalization of the correlation functions in a pro- 
cess in which the norm depletes. Of course, such a situation can not be considered as aging 



behavior. (In fact, in a recent paper |T^ this effect has even been proposed as a general 
explanation for aging.) 

III. CONCLUSIONS AND ACKNOWLEDGMENTS 

We have shown that a simple, solvable model of a random walk near a wall can exhibit 
aging behavior which illuminates many features that lead to aging behavior in more com- 



plicated (disordered) systems, using a domain-growth picture 0. Of course, systems with 
many interacting degrees of freedom such as spin glasses or folding proteins exhibit a non- 
trivial phase space structure [Jl2[ which leads to slow relaxation and aging behavior. Their 
dynamics is described merely on a coarse, phenomenological level by such a simple model. 
But the connection between the micro-dynamics and the macroscopic phenomena is not 
only beyond the scope of this article, but as well itself very much under development still. 
Instead of "explaining" experimental or theoretical results, this random walk model is meant 
to illustrate some of the questions involved. (After all, it is still rare to consider situations 
with broken time-translational invariance, and thus violations of the fluctuation dissipation 
theorem, without which two-time correlation functions would be redundant.) Furthermore, 
we have discussed some of the pitfalls in identifying aging behavior in systems which do not 
conserve the norm, and how to obtain the intrinsic features of such systems. 

1 would like to thank Maya Paczuski for discussing some of the issues considered, and 
Eli Ben-Naim for a critical reading of the manuscript. 
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FIGURES 
FIG. 1. Log-log plot of the two-time correlation function C{t + ty^^t^) (arbitrary scale) from 

numerical simulations of a walk near a reflecting wall. Each plot contains data for 2*~^ < t«, < 2* 

for 2 < j < 14 where i labels each graph from bottom to top. Initially, for t < tuj, each correlation 

function falls like l/\/t with a crossover at t ^ t^, after which the walk notices the effect of the 

(domain-)wall and the function falls like 2/\/t for t ^ t^. The continuous line is inserted to guide 

the eye and show that the plot of C for each t^, indeed falls like t^^''^ in two separate regimes which 

only differ by a factor of 2. 

FIG. 2. Log-log plot of the normalized, intrinsic two-time correlation function C"^*'"(t-|-t^, t^) 
from numerical simulations of a walk near an absorbing wall. (Each plot is shifted to avoid overlaps.) 
Each plot contains data for 8*^^ < i«; < 8* for 1 < i < 5 where i labels each graph from bottom to 
top. Initially, for t < tw, each correlation function falls like l/\/t with a crossover at t ^ t^, after 
which the effect of the (domain-)wall becomes noticeable and the function falls like 1/t. 

FIG. 3. Scaling plot f{t/tyj) ~ \/i C'°*'^(t-|-tu,, t^) as a function of the scaling variable t/tw for 
the data in Fig. 2. All data collapses reasonably well onto a single scaling graph which is constant 
for small argument and falls like an inverse square-root for large argument (such as the dashed line 
drawn for reference). 
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